One-dimensional infinite component vector spin glass with long-range interactions 
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We investigate zero and finite temperature properties of the one-dimensional spin-glass model for 
vector spins in the limit of an infinite number m of spin components where the interactions decay 
with a power, a, of the distance. A diluted version of this model is also studied, but found to deviate 
significantly from the fully connected model. At zero temperature, defect energies are determined 
from the difference in ground-state energies between systems with periodic and antiperiodic bound- 
ary conditions to determine the dependence of the defect-energy exponent 9 on a. A good fit to 
this dependence is 8 = | — a. This implies that the upper critical value of a is 3/4, corresponding 
to the lower critical dimension in the rf-dimensional short-range version of the model. For finite 
temperatures the large m saddle-point equations are solved self-consistently which gives access to 
the correlation function, the order parameter and the spin-glass susceptibility. Special attention 
is paid to the different forms of finite-size scaling effects below and above the lower critical value, 
a = 5/8, which corresponds to the upper critical dimension 8 of the hypercubic short-range model. 

PACS numbers: 75.50.Lk, 64.60. F-, 02.60.Pn 



I. INTRODUCTION 

The problem of understanding the physics of spin 
glasses in the form of simple model systems incorporating 
frustration and random disorder has challenged theoret- 
ical physicists for the last forty years 1 . Although signif- 
icant progress has been made 2 , predominantly through 
extensive numerical simulations, a number of important 
puzzles are still unsolved, and we do not have a clear 
understanding yet of the nature of the spin-glass phase 
and whether replica-symmetry breaking, the hallmark of 
the spin-glass state in mean- field theory 3 , carries over to 
systems in low dimensions d. A 

Unlike the case of ferromagnets, we have currently 
no means of performing a well-behaved perturbative ex- 
pansion of the replica-field theory of spin glasses in di- 
mensions below their upper critical dimension d u = 6. 
While this program effectively starts from the d = oo 
Sherrington-Kirkpatrick (SK) 5 model to understand be- 
havior in finite dimensions, an alternative approach is to 
consider vector spin glasses with an infinite number of 
spin components m — oo, but arbitrary spatial dimen- 
sion d, and extend these results in a 1/m expansion to 
the physically more relevant cases with finite to. 6 This 
approach appears particularly suitable as the to = oo 
limit implies a number of simplifications as compared to 
the models with finite to. Firstly, the model is replica 
symmetric even in the mean-field limit 7 , in contrast to 
the Ising, XY or Heisenberg spin glasses usually con- 
sidered (with m = 1,2 and 3 spin components, respec- 
tively). Secondly, it is tractable by analytical and nu- 
merical means. For calculations at zero temperature, it 
turns out to be very useful that the metastability afflict- 
ing finite-m spin glasses disappears 8 ' 9 , making it numer- 
ically straightforward to determine ground states. At 
finite temperatures, the to — > oo limit leads to saddle- 



point equations which allow for the exact calculation of 
correlation functions of finite samples for both mean-field 
and non-mean-field models 10 . 

On the other hand, the to = oo model has peculiari- 
ties. It has been shown that the upper critical dimension, 
d u , which is six for spin glasses with a finite number of 
spin components, is elevated to eight 11 . Likewise, the 
lower critical dimension at which a finite-temperature 
transition first occurs, appears to be also increased from 
that of systems with a finite number of spin components. 
It has been estimated from numerical studies that 12 ' 13 
di w 6 . The mechanism of effective dimensional reduc- 
tion that is at work in lifting d u also leads to a violation of 
hyperscaling even below the upper critical dimension 11 . 
Finally, regarding the numerical calculations considered 
here, one should note that an order of limits, m — > oo 
before N — > oo, is used which is opposite to that used 
in field theoretic calculations 11,14 . Taking the infinite- 
component limit first might be considered the zeroth- 
order term in a 1 /to expansion around the field-theoretic 
calculation 15 . 

From studies of ferromagnets it was realized many 
years ago 16 that systems with long-range, power-law in- 
teractions in low dimensions could be used to model the 
non-trivial critical behavior of the kind expected in short- 
range systems of higher dimensions. Similar observations 
were later made for spin glasses 17,18 . It was subsequently 
realized that such models are useful for numerical stud- 
ies, as finite-size corrections, known to be strong for spin- 
glass systems, depend on the linear extension of the lat- 
tice. Hence, studying, for instance, a one-dimensional 
system with interactions which fall off with distance with 
a power a allows one to access significantly larger (linear) 
system sizes than studying similar systems on hypercu- 
bic lattices 19-21 . While these first works considered Ising 
spin glasses, Potts 22 and Heisenberg 23-25 models have 
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FIG. 1. (Color online). The Id power-law spin-glass model 
on a ring geometry. The left panel shows the fully connected 
version where the magnitude of the interaction strength falls 
off with distance. The right panel shows the diluted model 
with the bond existence probability falling off with distance, 
whereas the bond strengths are distance independent. 

also recently received some attention. For the m = oo 
limit considered here, this approach appears to be well 
suited, as reasonable system sizes in the dimensions d > 6 
where a finite temperature phase transition occurs are 
nearly inaccessible with current computational resources. 
For the case of Ising spin glasses, diluted lattices have 
been used to reach even larger linear dimensions 26-28 , 
and the usefulness of this approach for the m — oo model 
will be discussed in some detail below. 

The paper is organized as follows. In Sec. II the 
model and some theoretical preliminaries will be intro- 
duced. Section III discusses the phase diagram of the 
one-dimensional spin glass with power-law interactions 
for finite and infinite m, and the scaling and finite-size 
scaling in the vicinity of the critical point. In Sec. IV, 
we report on the results of ground-state calculations and 
a study of the defect energies. The critical behavior is 
examined with finite-temperature methods in Sec. V. Fi- 
nally, Sec. VI contains our conclusions. 

II. THE MODEL 

In this paper we study flavors of the well-known 
Edwards- Anderson Hamiltonian 

1 L 

where the € K m , i = 1, . . . , L, are vector spins with 
m components, normalized as \Si\ = \pm. The spins are 
organized in an effectively one-dimensional (Id) geome- 
try, either chosen to be a chain with periodic boundary 
conditions or a ring as depicted in Fig. 1. In this fully 
connected version the exchange interactions are between 
all spin pairs, decaying as a power law with distance, 



where — | — | and Lpij is a standard normal random 
variable. 
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FIG. 2. (Color online). Correspondence between the Id spin- 
glass model with power-law interactions characterized by an 
exponent a and the short-range model on a hypercubic lat- 
tice of dimension d. The upper part of the figure applies to 
finite spin dimensions m, whereas the lower part describes 
the limit m — oo discussed here. Increasing a corresponds to 
decreasing the analogous lattice dimension d. 

For the case of Ising spins (m = i) 17 - 2 °. 26 > 27 > 29 an d, 
more recently, Heisenberg spins (m = 3) 24 ' 25 , this model 
has been extensively investigated. It is found that, as the 
range of interactions is tuned by varying <r, the model 
has a behavior which mimics that of the short-range spin 
glass as its dimension d is tuned. For large a the spin- 
glass transition temperature Tsg = 0. This corresponds 
to dimensions below the lower critical dimension di. A 
non-mean field regime is adjacent at intermediate a con- 
tinuing on to a mean-field region for small cr, which cor- 
responds to the dimensions above the upper critical di- 
mension, d u . Finally the SK model 5 is obviously reached 
in the limit a — > 0. The phase diagram of the model as a 
function of a is shown in Fig. 2, and will be discussed in 
more detail below. A dictionary can be set up relating 
the behavior of the Id power-law model at a given a and 
a corresponding short-range model on hypercubic lattices 
of dimension d. 26,27 We shall see here that for the limit 
of an infinite number m of spin components, the phase 
diagram is modified as shown in the lower part of Fig. 2. 

A. Choice of couplings 

We studied the model (1) in different variants and on 
different geometries. 

Fully connected model 

The fully connected system implied by Eq. (1) is real- 
ized with interaction constants 

J tJ=c (a,L)^, (3) 

where ipu G A/"(0, 1) are standard normal random vari- 
ables. It is the strength of interactions that falls off as 
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1 jr a here. The mean- field transition temperature, 
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diverges for a < 1/2, unless we prevent this by an ap- 
propriate L dependent choice of the normalization factor 
c(cr, L), for instance by requiring that 
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(c) 



1. 



(5) 



which fixes c(er, L). While this is only strictly necessary 
for a < 1/2, we apply the same normalization for all a. 
Clearly, the limit a — > corresponds to the SK model. In 
fact, it can be shown that mean- field theory is exact (at 
any temperature) for all a < 1/2. 30 For numerical sim- 
ulations employing single-spin manipulations, this fully 
connected model is slow as the number of bonds equals 
L(L— l)/2, so that the cost of a lattice sweep of updates 
scales quadratically with the system size L. 



Bond- diluted model 

To improve on this costly update for the fully con- 
nected model and allow numerical studies to get closer 
to the large system limit, a number of authors have con- 
sidered a diluted version of the Id power-law spin glass 26 . 
Its Hamiltonian reads 



1 



L 



(6) 



where now Jy € A/"(0, 1), but the probability distribution 
of the dilution variables g {0, 1} falls off with the 



distance as 




P <Pij, 
otherwise, 



(7) 
(8) 



with p £ W[0,1] a uniform random number from the in- 
terval [0, 1]. To ensure that the form of pij is a proper 
probability density function, we normalize 2 



Pij = 1 - exp(-A/r 2 /) 



(9) 



and determine A by fixing the average coordination num- 
ber 



L-l 
i=l 



xL- 



(10) 



Unless stated otherwise, for the data discussed here we 
used z — 12, corresponding to a hypercubic lattice at 
the probable lower critical dimension d[ — 6 13 . We ap- 
ply the Newton method 31 in A in order to iterate the 
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3. (Color online). The non-random part c(a, L) /rfj 
the interaction constants in Eq. (3), i.e. for the fully con- 
nected model, for the line and ring geometries and the re- 
summed (Hurwitz) version (a — 3/4 and L = 256). Note the 
logarithmic scale of the ordinate. 



probabilities PiL until their sum equals the desired coor- 
dination number to a certain precision. The factor l/\/~z 
in Eq. (6) ensures that TgQ F = 1, consistent with the 
fully connected model. This diluted version of the model 
was previously studied for the Ising 27 , Heisenberg 25 and 
p-spin 32 spin-glass models. The authors of Ref. 26 and 
subsequent studies claimed this model to be in the same 
universality class as the fully connected system. As we 
will see below, however, this is not the case for a > 1. 
Also, certain properties differ for a < 1/2. In the limit 
a — > 0, the diluted system corresponds to the Viana- 
Bray (VB) model 33 . Numerically, the diluted system 
with zL/2 bonds reduces the sweep time from quadratic 
to linear in L. 



B. Choice of geometry 

Two different effectively one-dimensional geometries 
have been previously considered in studying power-law 
spin glasses: a ring of spins 19 as depicted in Fig. 1 and the 
possibly more natural linear chain with periodic bound- 
ary conditions without any embedding space 34 . In the 
ring model, the distances are measured according to the 
Euclidean metric in the plane, 



L . ( ir\i- j\ 
— sin 

7T V L 



(11) 



and periodic boundaries are incorporated automatically. 
In the chain formulation, spins are located at integer po- 
sitions i on a straight line with distances 



mm(\i-j\,L 



(12) 



again assuming periodic boundary conditions. While one 
expects the specific form of the geometry to influence 
the finite-size behavior, in the limit of large distances on 
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large chains or rings, both formulations become equiva- 
lent. Universal properties, of course, should not depend 
on these details. On the other hand, one might argue 
that finite-size corrections, which are notoriously impor- 
tant in the study of spin-glass systems, will differ between 
the two formulations and might thus lead to an effective 
advantage for one or the other form. 

When studying long-range interactions one needs to 
be careful about defining a controlled approach to the 
thermodynamic limit. Using the periodic boundary con- 
ditions preferred to suppress boundary effects, each spin 
effectively interacts with an infinite set of periodic im- 
ages. The resulting infinite sums are usually performed 
in reciprocal space (Ewald summation). For the fully 
connected one-dimensional chain, i.e., Eq. (12), they can 
be performed without cut-off. Summing over images for 
the couplings of Eq. (3) one obtains the effective coupling 



Jij = c(a,L), 
c(a,L) 
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with the Hurwitz Zeta function 35 
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(14) 



The corresponding mean-field critical temperature is 
then 



c{a,Lf 



12,12,7+1 



E 



,(15) 



which, with the normalization TgQ F (c) = 1, fixes c(a,L). 

While for ferromagnetic systems summation over im- 
age "charges" is crucial, for a spin-glass system with av- 
erage magnetization (to) = it should not change the 
asymptotic behavior 36 , ft might modify finite-size cor- 
rections, however. In Fig. 3 we compare the non-random 
part of the interactions for the ring geometry as well as 
the bare and summed chain interactions. The differences 
away from r — are small and, as we shall see below, 
the alterations of the finite-size scaling (FSS) behavior 
are rather minor. Note that for the special case a = 1 
the constants c(a, L) for the ring and the summed line 
geometries coincide, which is easily understood from the 
identity 35 



C(2,ri i /L)+C(2,l-rl J ./L) 



Unless stated otherwise, all of the calculations presented 
below have been performed for the ring geometry. 



III. PHASE DIAGRAM AND CRITICAL 
BEHAVIOR 

To understand the expected critical behavior of the 
model in the to = oo limit, it is useful to review and 
generalize the results for the to < oo case. The most 
distinct feature of the to = oo limit on hypercubic lattices 
is the elevation of the upper critical dimension (UCD) to 
d u = 8 and the accompanying violation of hyperscaling 11 . 



A. Mean-field critical exponents 

Recall that the mean-field exponents of (say) the Ising 
ferromagnet are 

a = 0, p = 1/2, 7 = 1, v = l/2, 77 = 0. (16) 

These satisfy hyperscaling, 

dv = 2 - a, (17) 

(only) at the standard upper critical dimension d u = 4. 

For the (to < oo) spin glass, the upper critical dimen- 
sion is 1 d u — 6 and the mean-field value of the exponents 
are 

a = -l, 13 = 1, 7 = 1, v = \/2, r? = 0. (18) 

These again satisfy hyperscaling only at the upper critical 
dimension, d u = 6. a and /3 are the exponents of the SK 
model 1 . 

For the m = oo spin glass the upper critical dimension 
is 11 d u — 8. This model, however, violates hyperscaling 
also in the non-mean-field regime, and the hyperscaling 
relation is replaced by a dimensionally reduced version, 



{d-2)v = 2- a. 



(19) 



The mean-field exponents for the m = oo model are the 
same as those of the Ising spin glass i.e. 

a = -l, = 1, 7 = 1, ^ = 1/2, v = 0- (20) 



B. Finite-size scaling above the upper critical 
dimension 

To understand the behavior of the model in the mean- 
field regime, it is useful to recall the relevant form of 
scaling and FSS above the UCD. Below the UCD, finite- 
size corrections depend on the ratio of correlation lengths 
in the finite and infinite systems, £l/£oo ~ Lt v , where t = 
(T — T c )/T c is the reduced temperature. For a singular 
quantity A, we therefore expect the FSS form 37 



A ~ L K ' V A{L X / V t), d<d n 



(21) 
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where k is the critical exponent associated to A. For a 
dimcnsionless quantity such as the finite-size correlation 
length normalized by the system size, we expect 



£/L~ X{l}l v t), d<d u . 



(22) 



At and above the UCD, FSS should hold with mean-field 

exponents with the role of the correlation length £l ~ L 

taken on by some effective length 38 ' 39 Cz, ~ l ~ JV 1/d " = 
L d/d u guch that 



A~N K / d " v A(L d / d ""t), d>d u , 



(23) 



where v and k take on their mean-field values. Similarly, 
for the case of a dimcnsionless quantity, we arrive at 



£/ L d/d u „ X{L d / d " v t), d>d u . 



(24) 



We can therefore extend the hyperscaling law beyond its 
usual range of validity d < d u by replacing the correlation 
length exponent v with a renormalized value 



v, d < d u , 
d u v/d — d u /2d, d > d u , 



(25) 



since then dv' = 2 — a in all dimensions. As a conse- 
quence, at d = d u we find 



du' = d u /2 = 2 - a, 



(26) 



leading to a — for d u = 4 and a = — 1 for d u = 6. 

Comparing the critical scaling of the Landau- 
Ginzburg- Wilson (GLW) effective Hamiltonian for the 
one-dimensional long-range model and the short-range 
model in general dimensions d, one infers that close to 
the UCD one has 26 



d e g 



2(7 - 1 



(27) 



As we will see below, v — 1/(2<t — 1) = d e s/2 in the 
mean-field region, such that d/d u = d e s/d u = vjZ and 
djd u v = 1/3 and it follows from Eqs. (23) and (24) that 



Xsg ~ NWCttNV*), 

<ZEA~iV- 1/3 Q(w 1/3 ), 



L v l' i 



(28a) 
(28b) 

(28c) 



C. Id long-range spin glass 

The long-range Ising spin glass was discussed analyt- 
ically in Refs. 17, 18, 40, and 41. In Rcf. 42 it was 
proven rigorously that there is no phase transition for 
a > 1. Studying the effective GLW Hamiltonian in 
replica space 17 , it was inferred that there is a finite- 
temperature phase transition for 1/2 < a < 1, which 
is of mean-field type for 1/2 < a < 2/3 and of non- 
mean-field type for 2/3 < a < 1. Therefore, the lower 



critical ai — 2/3 corresponds to the upper critical d u — 6 
for systems on hypercubic lattices and similarly for the 
upper critical a u = 1 and the lower critical di . It is useful 
to set up a dictionary of correspondences between the Id 
long-range model and the short-range models on hyper- 
cubic lattices, cf. Fig. 2. Arguments were given in Larson 
et al. 32 that the effective dimensionality for 2/3 < a < 1 
was approximately given by 27 



d e fi = 



2 - y(d cS ) 
2a-l : 



(29) 



where r/(d e f{) is the exponent of the corresponding short- 
range model. The upper critical a u = 1 can be inferred 
from the result 18 
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LR 



1 - a 



(30) 



for the long-range defect-energy exponent (see Sec. IV B). 
Since 8sr = — 1 for the short-range Ising spin glass in Id, 
a finite-temperature transition ceases to exist at 6*lr = 0, 
that is, at a u = 1. 

In order to determine the correlation length £ for long- 
range models with power-law interactions, one uses the 
fact that the propagator is modified from the well-known 
Ornstein-Zernicke form to (in reciprocal space) 43 ' 44 



G(k) 



1 

to 2 + k 2a - 1 ' 



(31) 



For the spin-glass, the analogous form is for the spin-glass 
correlator xsc(fc)- Consequently the second- moment 
definition of the correlation length is modified to 



1 



2sin(fc mi „/2) 



XSG(O) 



- 1 



l/(2<x-l) 



(32) 



(2tt/L) Si e 



XSG(k n 

, where e\ is a lattice basis 



with k min 
vector. 

At criticality (where m = l/£ = 0), it is found that 
the Gaussian propagator <~ l/k 2a ^ 1 does not receive any 
corrections away from mean field 45 and hence 



2-77 = 2(7-1, (7 < 1. 



(33) 



The upper critical value o u = 1 also follows directly from 
observing that at the lower critical dimension (LCD), 
we expect the critical correlation function decay G(r) <~ 
l/r d_2+J) to be constant, i.e., di — 2 + r/ = 0. Since we 
have d = 1 and 2 — r\ = 2u — 1, it follows that a u = 1. 

In the mean-field regime of the Ising model, Kotliar, 
Anderson and Stein (KAS) find 17 

v = l/(2cr- 1), 1/2 < cr < 2/3. (34) 

This implies that 

7 = (2 - rfiv = 1, 1/2 < cr < 2/3. (35) 

Using modified hyperscaling, Eq. (26), for cr < 2/3, we 
expect 



dv' = 1 • 



6 



2cr-l 2/(2(7-1) 



3 = 2- a, 



(36) 
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i.e., a = — 1 and v' = 3 for 1/2 < cr < 2/3 and, conse- 
quently, 

^ = (2 - a) - 7 = ^ 1/2 < (T < 2/3 ^ (37) 

Note that with 

2-tj = 2tr-l = (2- r? MF )Kff (38) 

and 

1 



2cr- 1 



^MF^cff, 



(39) 



where (i e ff = 2/(2cr — 1), all exponents take their expected 
mean-field values. 

In the non-mean- field regime, 2/3 < a < 1, KAS 
showed that an expansion around the lower critical value 
er; = 2/3 in the variable e = a — 2/3 was possible and 
yielded to first order in e 



v 3 



(40) 



Expansions around the upper critical a u = 1 have also 
been proposed 17 ' 41 . 



D. Id long-range m — oo spin glass 

For the m = oo model we know 11 that the UCD is 
elevated from the usual (spin-glass) d u = 6 to d u = 8 
and that, additionally, there is a failure of hyperscaling, 
even below the UCD. For hypercubic lattices at the UCD, 
FSS should work in N (see above), e.g., 



Xsg 



L1 / v = L 2- n = N (2-r,)/d = N i/i^ d > 8 ( 41 ) 



For the Id long-range model, we expect the long-range 
form of the exponent of the correlation function, Eq. (38), 
to carry over to the m — oo model. At the lower critical 
(7j , where mean- field behavior first becomes modified, we 
should find 



(2 - ri)/d = 2 - r] = 2a t - 1 = 1/4 



(42) 



or <t; = 5/8. Therefore, the mean-field regime is here 
defined as 



1/2 < er < 5/8. 



(43) 



Of course, this range can be also obtained directly via 
the calculational methods in Green et a/. 11 . The effec- 
tive correlation length exponent thus becomes v' = 4 
and, due to dimensional reduction, we expect a modified 
hyperscaling relation to hold, 



(d-ey = 2- 



(44) 



with some violation-of-hyperscaling exponent for the 
long-range case. Since we should have a = — 1, we infer 



9 = 1/4 for 1/2 < a < 5/8. The "bare" correlation 
length exponent should be unaltered, 



2(7 - 1 



(45) 



so that then 7=1 and (3 = 1 as expected from mean-field 
theory. 

The FSS forms of the critical quantities become modi- 
fied by d u — 8 according to the discussion outlined above 
to read 



<Zea ~ N-^QitN 1 / 4 ), 



XitN 1 / 4 ). 



(46a) 
(46b) 

(46c) 



A consistent definition of the violation-of-hyperscaling 
exponent is given by 
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2a-l = 2/d eS , 5/8 < cr, 

1/4, 1/2 < a < 5/8. 



(47) 



When a > 5/8, this follows from the form of the prop- 
agators at T c , which go as l/fc 2cr_1 , and the results in 
Rcf. 11. is used in scaling relations which involve the 
dimensionality d when one replaces d by d— 0. Thus the 
scaling relation f3jv = (d — 2 + r])/2, with the replace- 
ments d ^ d — 0, 2 — r] = 2a — 1, and d = 1 becomes 



P/v= (3-4ct)/2. 



(48) 



This is consistent with our numerical results shown in 
Fig. 25. It is possible to determine exactly the value of 
the upper critical a u from generalizing the argument that 
at the lower critical dimension d — 2 + r) = 0. Replacing 
once again d — > d — and 2 — rj = 2a — 1 , one has at 
a = a u , 1 - (2a u - 1) - (2a u - 1) = 3 - Aa u = 0, so 

o u = 3/4. 

As will be discussed below, the defect-energy calcu- 
lations for m = oo presented here, cf. Fig. 12, can be 
summarized as 9(a — 3/4) = 0, i.e., a u — 3/4, and 
9(a = 1/2) = 1/4, which lead us to conjecture that 



'LR 



3 
4 



(49) 



In the following, wc refer to #lr simply as 9. The form 
(49) works over a rather wide range of cr, even for cr > 3/4 
in the fully connected model. When a > 3/4, there is no 
finite temperature transition. We have been unable to 
give a formal derivation of this result, but suspect that 
this might be possible by generalizing the formalism of 
Aspelmeier et al. 46 to spatially varying solutions. One 
can however understand Eq. (49) from the following con- 
siderations. From the scaling arguments in Ref. 18 for 
long-range Ising spin glasses 29 = 2d — 2a. This would 
also follow from the formalism of Ref. 46 which would re- 
sult in an expression for the variance of the defect energy, 
(which scales as L 26 ) proportional to a double sum over i 
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and j of [J^Jav if a spatially non-uniform solution for the 
defect energies is studied. We have to consider how the 
failure of hyperscaling for the large to limit might affect 
this relation. 

is not a critical point exponent, but an exponent asso- 
ciated with the fixed point at zero temperature. For it, wc 
suspect that the mean-field form of 6 = 1/4 is relevant 
for both a greater than and less than 5/8, since zero- 
temperature exponents like can usually be obtained 
by a simple minimization of the defect energy, just as 
one determines mean-field behavior by minimizing the 
total energy of the system. Thus allowing for the fail- 
ure of hyperscaling, the equation 20 — 2d — 2a becomes 
26 = 2{d - 6) - 2a. With 6 = 1/4, Eq. (49) for 9 is ob- 
tained on setting d = 1. We would expect this argument 
to still be valid in the fully connected model even for 
a > 3/4 when there is no finite temperature transition. 

McMillan 40 ' 47 has argued that the relevant renormal- 
ization group equation for the flow of the temperature T 
near the lower critical dimension is 



dT 
dlni 



= -6T + cT 3 + .... 



(50) 



For 9 small and positive, (i.e., for a below, but close to 
3/4), one finds a fixed point at Tsg °c V9 cx \/3 — 4a. If 
we choose the proportionality constant so that Tsg(ct < 
1/2) = 1, then 



Tsg = V3 - 4a, 



(51) 



which fits the critical temperature quite well in the whole 
regime 1/2 < a < 3/4 (cf. Fig. 25). The eigenvalue at 
this fixed point is 



1 



26» 3 - 4a 



(52) 



This appears to be consistent with the data shown in 
Fig. 23 for the regime 5/8 < a < 3/4, but it cannot be 
regarded as anything but an interpolation formula, exact 
only at the end points a — 5/8 and a = 3/4. For a > 3/4 



48 „ = 

From the scaling relation 7 = v(2 — r\) with v given by 



T SG = 0. Then one expects 4 * z/ = -1/0. 



Eq. (52) one has the approximate result that 



4cr-2 



7 = 7 



3 -4a' 



(53) 



in the regime 5/8<a<3/4. 7=1 is thus expected at 
a = 5/8 and 7 — » 00 as a — > 3/4 from below. 

By combining Eq. (48) with Eq. (52) one finds j3 — 
1 throughout the interval 5/8 < a < 3/4. Thus the 
expectation is that (3 remains close to its mean-field value 
even in the non-mean-field region. 



IV. ZERO-TEMPERATURE CALCULATIONS 

We start our numerical investigations of the m = 00 
spin glass by studying its ground-state properties as a 



function of a. Due to the possibility of studying large 
system sizes, we first concentrated on the diluted model 
of Eq. (6). For most calculations, the ring arrangement 
was used. Compared to the ground-state problem for 
generic spin glasses which is found to be NP hard 49 , 
ground states for the m — 00 limit are much easier to de- 
termine. Starting out with the Ising model with to = 1, 
with increasing spin dimension the energy landscape sim- 
plifies gradually until, for to — > 00, all mctastability has 
vanished and the ground state becomes unique. This fact 
was already exploited for hypercubic systems in Ref. 13, 
where the lower critical dimension was determined with 
the defect-energy approach and for different boundary 
conditions. 

The possibility to realize the limit to — > 00 in numerical 
calculations rests on the fact that for a finite system of 
TV spins, the ground state occupies a finite dimensional 
sub-manifold in spin space 50 ' 51 , the dimension of which 
is limited by the rigorous upper bound 



{N)= (VSN + 1 - l) /2 ~ JV", n= 1/2, (54) 



where \x\ stands for the largest integer smaller than or 
equal to x. Hence, for each system size a finite num- 
ber m*(N) < m laax (N) of spin components is suffi- 
cient to describe the to = 00 model. For commonly 
used spin-glass models, the scaling is in fact weaker than 
m*(N) <~ N 1/2 . For the SK model realized, e.g., in the 
limit a — > of our fully connected Id spin glass, one 
finds 6,50 fi — 2/5. As the degree of connectivity is low- 
ered, (i is reduced. We determined the required num- 
ber of spin components m*(N) for each single realiza- 
tion of the bonds Jij and computed the disorder average 
mo = [m*(JV)] av . The values of m*(N) are found to vary 
only slightly between disorder realizations, such that us- 
ing TO ac t ~ 1.1 [to* (iV)] av was sufficient to ensure that 
to = 00 ground states are found for all realizations. The 
procedure of determining the number of necessary spin 
dimensions to* will be described at the end of Sec. IV A. 

Due to the lack of metastability for m — > 00, it is quite 
straightforward to determine true ground states numer- 
ically. Here, we employ a local spin-quench procedure, 
for which the spins are iteratively aligned with their re- 
spective local molecular fields H i7 so that the new value 
of the spin Sj is given by 



s: 



J€JV(») 



Jij Sj , 



(55) 



where the sum runs over the set Af(i) of connected neigh- 
bors of the spin at site i. It is easily seen that alignment 
of each spin with its molecular field is a necessary con- 
dition for the system to be in its ground state. For the 
present case of a system without metastable states 50 , it 
is also sufficient. These updates are interspersed with 
sweeps of over-relaxation moves to speed up convergence, 
which have also been found to improve the decorrelation 
of systems with finite spin dimension m in Monte Carlo 



simulations 52 . These moves, again being local, preserve 
the energy of the whole spin configuration since the up- 
dated spin is merely rotated around its local field and 
therefore moves at constant energy. The simplest way 
of implementing such a procedure, in particular for the 
case of arbitrary spin dimensions m, is to reflect the spin 
along Hi, such that 



S' 



, Si • Hi 



H, 



(56) 



This maximal movement can also be argued to lead to a 
maximal decorrelation effect within the constant-energy 
manifold of single-spin movements. The whole procedure 
of spin-quench and over-relaxation moves can be imple- 
mented very efficiently, since only a few elementary oper- 
ations are required for each step, and no random numbers 
are involved. 



A. Ground-state properties 

Finite-size corrections to the ground-state energy 
of spin glasses have been extensively discussed re- 
cently for the case of the short-range Edwards- Anderson 
system 53-55 , spin glasses on the Bethe lattice and random 
graphs 56 , and the SK model 57 . The dominant contribu- 
tion for short-range systems is due to the presence of 
domain-wall defects, leading to corrections proportional 
to 54 L d ~ e . For the system studied here, however, these 
effects, although presumably present, are masked by cor- 
rections stemming from the power-law nature of the in- 
teractions. As indicated in Eq. (4), the relevant energy 
scale for the case of unrenormalized coupling strengths, 
i. e., c(ct, L) = 1, is set by the integral over the couplings, 

L-l 



r 

i=i » 



Hence, the ground-state energy per spin and spin- 
component, i.e., e{L) = E/Lm, is expected to scale as 



e'(L) 



c'L 1 



(57) 



where further finite-size corrections stemming from the 
presence of domain-wall excitations etc. have been ne- 
glected. We will see below that these are sub-leading 
and cannot be resolved by the numerics. The primed 
quantities in Eq. (57) are meant to indicate the unrenor- 
malized case with c(a,L) = 1. For large systems, we 
therefore expect different limiting behaviors depending 
on whether a ^ 1/2, viz. 



e'(L) 



2e^ 



- 1-2(T 



LV2- 



cr > 1/2, L > 1, 

a < 1/2, L > 1, 
(58) 
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FIG. 4. (Color online). Average ground-state energies of the 
Id power-law model in the m — ¥ oo limit. Top panel: ground- 
state energies e as a function of system size L and interaction 
range a. The lines show fits of the form (61) to the data. 
The corresponding correction exponents b are shown in the 
middle panel. The blue lines indicate the expectations for 
the unrenormalized energies e'(L) which are b' — a — 1/2 for 
a < 1/2 and b' = 2a — 1 for a > 1/2. The bottom panel 
shows the resulting asymptotic ground-state energies eoo, for 
the diluted and the fully connected model, respectively. 
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with logarithmic scaling right at a = 1/2. If we choose to 
make the energy scale convergent for a < 1/2 by setting 



L-l 



^) 2 = Ei. 



V7 
= 1 *J 



(59) 



we instead consider e(L) — e'(L) ■ c(a,L) with limiting 
behavior 



e(L) 



cL 1 ~ 2a , a > 1/2, L> 1, 
const, a < 1/2, L > 1. 



(60) 



For the case of the diluted model, similar considerations 
lead to the same results, where now A of Eq. (9) takes 
on the role of c(a,L) 2 . In Fig. 4 we present the results 
of the scaling of ground-state energies. The correction 
exponents result from fits of the general form 



- -b 



e(L) — + cU 
e'(L) = e' 00 + c'L- b ' 1 



(61) 



to the data. The number of disorder realizations used for 
the ground-state calculations are summarized in Table I. 
As is seen from the middle panel of Fig. 4, the predictions 
b' = 2a - 1 for a > 1/2 and V = a - 1/2 for a < 1/2 for 
the e'(L) and b — 2a — 1 for a > 1/2 in the renormal- 
ized case are borne out well in the data. For a < 1/2, 
where we predict b — 0, sub- leading corrections become 
visible. The resulting correction exponent b = 2/5 is 
consistent with the expectations for the SK model, cf. 
Ref. 58 and the discussion in Sec. VD below. For the 
renormalized energies e(L), we see a dip of the correc- 
tion exponent for 0.5 < a < 0.6, which is possibly due to 
additional finite-size effects resulting from the crossover 
between the forms for a ^ 1/2. As shown in the bottom 
panel of Fig. 4, the asymptotic ground-state energy eoo 
smoothly increases for interaction ranges a > 1/2. For 
a < 1/2 it is independent of a and takes the value —1 in 
the fully connected version of the model 8 . The indepen- 
dence of this non-universal quantity on a in this regime 
is a clear sign of the exactness of mean-field theory for 
a < 1/2 as proposed in Ref. 30. For models of lower 
connectivity, however, this energy is increased. Calcula- 
tions on a Bethe lattice 58 are consistent with our results, 
for instance e oa (a = 0.1) = —0.8784(1), if the average 
coordination number is taken into account. 

At T = the free energy F reduces to the internal 
energy. We can use this fact to consider its sample-to- 
sample fluctuations, 



(62) 
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FIG. 5. The scaling behavior of the sample-to-sample fluctu- 
ations characterized by the exponent Qf for a = 0.1. This 
showcase example for the fully connected model indicates 
the correctness of the prediction Of = 1/5 in the SK limit 
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FIG. 6. (Color online). Sample-to-sample fluctuations differ 
for the diluted and the fully connected model. Whereas the 
characteristic exponent is fixed to 1/2 for the diluted model, 
it has a jump from 1/5 to 1/2 when changing from a < 1/2 
to a > 1/2 for the fully connected model. 



in a FSS analysis. The scaling of this quantity has 
been the subject of a number of recent analytical and 
numerical studies see, e.g., Refs. 53, 59-63. For the 
(m = 1) Ising SK model, there has been some debate 
as to whether 64 6/ = 1/6 or 46 6/ = 1/4, but now 
there is growing consensus that 57 8/ = 1/6. For the 
7tl — oo model discussed here, the situation is less well 



studied. For the replica symmetric spherical SK model 
it was shown 59 that 8* = 1/3. As the m — oo limit is 
also replica symmetric 7 , one might suspect this value to 
carry over to the present case. In Ref. 63, however, it 
was argued on the basis of connections of the problem to 
bond chaos that, instead, Of = 1/5. For the long-range 
model studied here, we therefore expect 0/ = 1/5 in the 
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infinite-range regime er < 1/2 and a trivial 0/ = 1/2 for 
a > 1/2. Using our data for the ground-state energies 
to study this problem, we find clear power-law scaling 
of the distribution widths o~n irrespective of the chosen 
value of a. An example is presented in Fig. 5 for a = 0.1. 
As shown in Fig. 6, the value of 0/ for the fully con- 
nected model is found to approach Of = 1/5 in the SK 
limit as predicted, hence confirming that the m = oo 
model is not in the same universality class as the spheri- 
cal spin-glass. In the short-range limit a — > oo, we arrive 
at a trivial Of = 1/2. The transition between these 
two extremes, expected to be sharp at a = 1/2 in the 
thermodynamic limit, is found to be rather smeared out, 
however, indicating the presence of strong corrections to 
scaling. In contrast, the situation for the diluted model 
is found to be rather different with the fluctuation expo- 
nent being compatible with Of = 1/2 irrespective of cr, 
cf. Fig. 6. This is some first evidence of a lack of univer- 
sality between the fully connected and diluted versions 
of the Id power-law model. In fact, a similar behavior 
was predicted in Ref. 62 for the case of diluted Ising sys- 
tems, where it was attributed to the local heterogeneities 
caused by the fluctuating coordination number. 

We also checked the distribution functions of the 
ground-state energies for the diluted as well as the fully 
connected model and each power-law exponent a. While 
these show a non-trivial form for the (Ising) SK model 60 , 
Gaussian distributions have been reported for short- 
range models 53 . For the Id Ising power-law chain, a 
crossover from Gaussian to non-trivial has been found 
on moving into the mean-field regime 20 . The distribu- 
tion of ground-state energies for the Id to = oo model is 
analyzed in Fig. 7. For all values of the power-law expo- 
nent a considered here (0.1 < a < 1) the distributions 
seem to be compatible with a Gaussian. This is indi- 
cated by a showcase example for cr = 0.1 and N = 32 678 
with a quantile-quantile plot as well as an analysis of the 
skewness of the distribution, estimated by 



1 Ns 



Ei - E 



(63) 



in dependence of the system size for several values of 
a spanning the mean-field as well as non-mean-ficld 
regimes. Here, N s denotes the number of disorder realiza- 



tions. E — (1/N S ) is the usual estimator for the 

expectation value (E) and a\ = (1/(N S — 1)) Ylf=i(Ej — 

E) 2 estimates the variance a 2 E . Comparing to the results 
of Ref. 20, it is worthwhile to note that lattice sizes used 
there were considerably smaller (L < 192) due to the fact 
that a Monte Carlo simulation was employed on the fully 
connected Id power-law model. As shown in the lower 
panel of Fig. 7 for such small systems we see a signifi- 
cant decay for increasing lattice size similar to the one 
reported in Ref. 20. The data shown stem from calcu- 
lations using the diluted version of the model. We also 
checked the fully connected version, however, and found 
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FIG. 7. (Color online). Evidence of the Gaussian nature of 
the ground-state energy distribution for a range of different 
values of a. The top panel shows a normal quantile-quantilc- 
plot in a showcase example for a — 0.1 and the largest avail- 
able lattice size N = 32 768. Sample quantiles do not deviate 
from the theoretical straight line assuming a Gaussian dis- 
tribution function. The lower panel shows the scaling of the 
(modulus of the) skewness according to Eq. (63) for a — 0.1, 
0.5, 0.75, and 1. It vanishes for all a verifying the Gaussian. 



identical results (not shown). 

For the case of the to = oo model considered here 
another zero-temperature property concerns the average 
number to = [ro*] av of spin components required to form 
the ground state. We determine to* for each realization 
by ordering all spin vectors of the ground state config- 
uration into an to x TV-matrix M — {Si, . . . , Sjv} and 
performing a singular value decomposition (SVD) in or- 
der to calculate the number of non-zero singular values, 
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FIG. 8. (Color online). Top panel: singular values of the spin 
matrix M for a sample with N = 4096, <j = 0.1 for a quench 
into the ground state. The "precision" refers to the relative 
change of the configurational energy after a fixed number of 
sweeps. Middle panel: average number of ground-state spin 
dimensions mo for a = 0.1. The line shows a fit of the form 
(64) to the data. Bottom panel: comparison of mo for the 
ring geometry and the line geometry with bare and resummcd 
interactions (Hurwitz) for a = 3/4. 



which is equal to the rank of the matrix. Since the rank 
of a matrix determines the number of linearly indepen- 
dent columns, i.e., spin- vectors, it is the desired num- 
ber m* . In practice, we monitor the size of all singular 
values of M online as the quench into the ground state 
proceeds. As the precision of the ground-state determi- 
nation is increased, those singular values that vanish in 
the exact ground state will scale to zero, whereas all the 
other singular values reach finite limiting values. This is 
illustrated for a sample of size N — 4096 and a = 0.1 in 
the top panel of Fig. 8. 

For the average number of spin components in the 
ground state, we assume the scaling form 



m 



const, x 



(64) 



where the additive correction c can account for the fact 
that for small systems mo will not scale to zero, but will 
be mo = 1 for tiny systems with N — 1 and N = 2. 
We present the results of this analysis for the infinite- 
range value a = 0.1 and the diluted model in the middle 
panel of Fig. 8. Since mean-field theory is exact there, 
we expect to see the value /i — 2/5 found for the SK 
model 6,50 , which is indeed borne out rather well. For 
a > 1/2, n continuously decreases below /i = 2/5, cf. the 
summary of our data for the diluted model in Fig. 9. A 
fit to a parabola yields /Lt(<r) = 0.3995 - 0.55(cr - 0.504) 2 . 
Note that since mo ~ and due to the boundedness 
of A* < 2 / 5 < 1 , ii; follows that m /N ~ N^ 1 -> in 
the thermodynamic limit. Hence, as T — > the spins 
condense into a subspace of vanishing relative size, just 
as in the more familiar Bose-Einstein condensation 6 . 

To check for the influence of the different geometries 
introduced above in Sec. II B on the scaling results, we 
performed some calculations for the bare and resummcd 
line geometries. The effect of these changes on m is il- 
lustrated for a — 3/4, where we expect the largest devia- 
tions, in the bottom panel of Fig. 8. We find a small 
overall shift in mo but, as expected, no change in fj,. 
Moreover, there appear to be no significant alterations 
of corrections to the leading FSS behavior. 



B. Defect energies 

The defect-energy approach 65 is widely used in study- 
ing systems with spin-glass phases. It is based on the as- 
sumption that the cost -E^ef of the insertion of a system- 
size defect into a state of the ordered phase scales as 66 



E dc { oc L 



(65) 



where 9 is known as the spin-stiffness exponent. Gener- 
alizing Peierls' argument for the stability of the ordered 
phase of a ferromagnet, one predicts Tgc = whenever 
9 < 0, whereas the ordered phase is stable at finite tem- 
peratures for 9 > 0. The limiting case 9 — corresponds 
to the LCD of the system. Additionally, for the case 
9 < with a zero-temperature transition, 9 is related to 



the correlation length exponent as 1 



()<> 



-1/9. 
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FIG. 9. (Color online). Variation of the exponent /j, with the 
interaction range a. The points are results of fits of the form 
Eq. (64) to the data for the average number of occupied spin 
components in the ground state for the diluted version of the 
model. Gray full points show the quality-of-fit parameter Q 
(right scale). The purple line corresponds to the functional 
form (i{a) = 0.3995 - 0.55(cr - 0.504) 2 . 

Numerically, defect energies are conventionally deter- 
mined by comparing ground states of systems with a 
pair of different boundary conditions (BCs) chosen such 
that the respective ground states must differ by a relative 
domain-wall type excitation. Then the defect energy cor- 
responds to the energy difference. The most commonly 
used set of such BCs are periodic and antiperiodic bound- 
aries. The defect energy of a given realization is then 

AE=\E AP -E P \, (66) 

where the modulus is required since, for symmetric cou- 
pling distributions, the two boundary conditions are sta- 
tistically equivalent. For the case of the long-range ring 
geometry considered here, a ground-state search is per- 
formed for the original coupling configuration, yielding 
Ep. In the second step, the boundary exchange couplings 
are flipped to the antiperiodic state by choosing one ar- 
bitrary nearest-neighbor pair (S a , S a +i), a € 1,...,N 
(without having them necessarily interact in the diluted 
version of the model) and changing the sign of all in- 
teraction constants Jij between spin Si and Sj for all 
i =/= j if the shorter path between those two spins falls 
on top of the path between S a and S a +i. 19 A second 
ground-state search for this altered configuration then 
yields E A p. There has been some discussion in the past 
about whether this setup is suitable for the case of contin- 
uous spins, since both periodic and antiperiodic bound- 
aries induce some defects, such that the energy differ- 
ence E P — Exp does not directly correspond to a de- 
fect energy 13 ' 67-69 . Therefore, boundary conditions that 



FIG. 10. (Color online). Convergence of the defect energy 
for the diluted model as a function of the precision of the 
ground-state energy calculation. For large a this should be 
monitored accurately. 

allow directly to measure the energy of a single defect 
seem preferable. For the case of the m = oo spin glass 
on hypercubic lattices we indeed found different stiffness 
exponents from such "domain wall" BCs 13 . For the Id 
long-range system, however, it is not obvious how to im- 
plement such alternative prescriptions. 

The resulting energy differences are averaged, 

E dci = [AEl v = l\E AP ~Ep\l v (67) 

to yield an estimate E de { of the defect energy. Ground 
states were computed for a in the interval 0.1 < a < 1 
for the diluted version of the model. For most of the val- 
ues of a a range of system sizes 512 < L < 32 768 has 
been considered with about 3 000 disorder realizations for 
larger and up to 10 000 samples for smaller system sizes, 
see the parameters collected in Table I. We ensured con- 
vergence by monitoring E do { as the ground-state quench 
proceeds, cf. Fig. 10. This is of particular importance for 
large cr, where choosing a fixed precision fails to produce 
converged results for sufficiently large systems. 

The averaged defect energies are shown for the avail- 
able interaction ranges a in Fig. 11. To extract the stiff- 
ness exponents, we performed fits of the functional form 

E dcf = aL (l + b/L) (68) 

to the data. Similarly to our experience from the hyper- 
cubic systems, we found this form to describe the cor- 
rections rather well 13 . While this form parametrizes the 
leading analytical correction for 9 < 0, a constant would 
be asymptotically dominant over b/L for 9 > 0. Fits in- 
cluding a constant but no 1/L correction, however, are 
not found to describe the data well for 9 > 0, such that 



13 




FIG. 11. (Color online). Defect energies for the diluted Id 
power-law m — oo spin glass as a function of lattice size L for a 
number of different interaction ranges a. The solid lines show 
fits of the functional form (68) to the data. The corresponding 
fit parameters are summarized in Table I. The defect-energy 
exponent 9 changes sign in the vicinity of a = 0.77, cf. Fig. 12. 



we stick with the form (68) for all ct. Any non-analytic 
corrections, if present, appear to be sub-leading. These 
fits are shown in Fig. 11 and the corresponding fit pa- 
rameters are collected in Table I. Fit qualities are found 
to be high throughout, indicating the suitability of the 
form chosen in Eq. (68). 

Figure 12 summarizes our results for 9 as a function 



TABLE I. Estimates of the spin stiffness exponent 9 in the 
diluted model resulting from fits of the functional form (68) 
to the data, including system sizes A^min < N < iV max . Q 
denotes the quality-of-fit. 
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of ct. The stiffness exponent 9 clearly becomes constant 
at a value compatible with 9 = 1/4 in the infinite-range 
regime a < 1/2. To determine the upper critical ct u 
where 9{a u ) = 0, we performed a linear fit to the results 
in the range 0.5 < ct < 0.9, resulting in an intercept 
of a u — 0.76(3). In connection with the observation of 
a linear behavior in the regime 1/2 < ct < 3/4, this is 
compatible with the conjectured form #lr. = 3/4 — ct, 
cf. Eq. (49). This behavior is clearly different from the 
corresponding Ising spin-glass model with 0lr = 1 — ct. 
Notice that the change to the value 1/4 at ct = 1/2 is not 
due to a failure of Eq. (49) for a < 1/2, but because for 
ct < 1/2 we have rescaled the bonds down by a factor of 
L 1 / 2 ^" . If we had not done that 9 would have continued 
to be fitted by Eq. (49). 

For larger ct, however, we observe clear deviations of 
the data for the diluted model from 9 = 3/4 — ct. In fact, 
the data for E^ef at ct = 1.3, not shown in Fig. 11 (but 
see the lower left panel of Fig. 13), show a strong down- 
ward curvature, more resembling an exponential decay. 
A closer look reveals that there is no universality between 
the diluted and fully connected model for a > 1, where 
the properties of the diluted graphs change significantly. 
As has been shown in Refs. 70 and 71, Id graphs defined 
by Eq. (9) always percolate for ct < 1/2 and they per- 
colate for sufficiently large A (namely, for any z > 1/2) 
in the regime 1/2 < ct < 1. In contrast, percolation is 
(asymptotically) absent for ct > 1. For the defect-energy 
calculations considered here, this means that such non- 
percolating samples contribute AE = to E^ c f, leading 
to much smaller averages than expected from the scaling 
Eq. (65). For the average coordination numbers con- 
sidered here, such breakdown of percolation is only ob- 
served for very large systems, mostly beyond the reach 
of our numerical calculations. If we remove the links 
up to a finite range, however, for instance all nearest- 
neighbor links, the remaining graph does not percolate 



TABLE II. Estimates of the spin stiffness exponent 9 in the 
fully connected model resulting from fits of the functional form 
(69) (for a < 0.75) and (68) (for a > 0.75), respectively, to 



the data. 
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FIG. 12. (Color online). The stiffness exponent 9 extracted 
from the diluted model as a function of the interaction range 
exponent a. See Table I for the corresponding fit parameters. 
9 changes sign for the critical value o u — 0.76(3) marked by 
the shaded area. Below a — 1/2, the spin stiffness exponent 
levels off at the infinite-range value 9 — 1/4. The dashed line 
denotes the conjectured form — 3/4 — a, cf. Eq. (49). The 
inset shows analogous results for the fully connected model 
and a wider range of a. 



parameters and results for this model are summarized in 
Table II. In contrast to the diluted model, a constant 
was found to be a good description of the leading scaling 
corrections for 9 > 0, such that we used the form 

E dei = aL e + c (69) 

for a < 3/4 and the form (68) for a > 3/4. Apart from 
the fact that only the fully connected model does repre- 
sent the long-range universality class for a > 1, we also 
find scaling corrections for a < 1 to be less pronounced 
there than for the diluted system, such that it appears 
questionable whether considering the diluted model of- 
fers a significant advantage in terms of the precision and 
accuracy of the final results. 

Note that the results for 9 as a is increased are in 
contrast to those for the Ising case, where 6>sr = — 1, 
so that the (fully connected) long-range Ising system is 
governed by short-range behavior for 6*sr > 0lr viz. a > 
2 (whereas the diluted Ising system would be truly short 
ranged already for a > 1). The m = oo model, instead, is 
truly long-range everywhere, and crossover to the 9sr = 
— 1 of the nearest-neighbor Id chain system is not seen. 
Instead, the limit a —¥ oo of our fully connected model 
corresponds to the Id ladder system with 9 9 = — oo. 



V. FINITE-TEMPERATURE CALCULATIONS 



for a > 1 already for moderate sizes such that the long- 
range nature is lost. An alternative way of understand- 
ing this phenomenon is to note that the diameters of the 
graphs considered here grow proportional to (\ogL) s for 
1/2 < a < 1, corresponding to an infinite-dimensional 
or small-world graph, whereas they grow proportional to 
L for a > 1, corresponding to a truly one-dimensional 
graph 72 . This explains the strong downwards deviations 
of 9 from the form 9 = 3/4 — a seen in Fig. 12 for a > 1. 
Right at a = 1, we expect non- universality with 9 de- 
pending on the average coordination number z which, in 
turn, is a function of the parameter A in Eq. (9). 

This is illustrated in Fig. 13 showing our estimates for 
8 from diluted lattices with different average coordination 
numbers z. While, for a < 1, the measured 9 is essen- 
tially independent of expected, there is a dramatic 
z dependence right at a = 1, with the non-universal 9(z) 
approaching the value expected for the fully connected 
model in the limit z — > oo. Similar z dependent critical 
behavior at a = 1 was also recently found for random 
walks on such graphs 73 . For a > 1, on the other hand, 
the curves -Edof cease to follow power laws and, instead, 
cross over to the 9 = — oo or exponential decay expected 
for short-range one-dimensional systems 9 . 

These deviations are specific to the diluted model: in- 
dependent ground-state calculations for the fully con- 
nected model, summarized in the inset of Fig. 12, are 
consistent with 9 = 3/4 — a also for cr > 1. The relevant 



The critical behavior which arises when a < 3/4 can 
only be studied with techniques appropriate to finite 
temperatures. While the ground-state calculations were 
greatly facilitated by the disappearance of metastability 
in the m — ¥ oo limit, the use of a saddle-point proce- 
dure permits the exact solution at finite temperatures. 
As will be discussed in this section, this approach leads 
to an iterative set of matrix equations that allow one to 
determine the thermally averaged spin-spin correlation 
function which, in turn, gives access to the Edwards- 
Anderson order parameter, the spin-glass susceptibility 
and the spin-glass correlation length. Due to the nature 
of the matrix equations, use of a diluted model does not 
have any computational advantages, so that all calcula- 
tions have been performed on the fully connected model 
of Eq. (3). 



A. Saddle-point equations 

The saddle-point equations for the m = oo model 
were first derived in Ref. 10 and later discussed again 
in Refs. 6, 9, and 50. Starting point is the partition func- 
tion corresponding to the Hamiltonian in Eq. (1), 
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FIG. 13. (Color online). Dependence of defect energies and 
stiffness exponents in the diluted model on the average coordi- 
nation number z. While for a < 1, this is merely a finite-size 
correction, we observe a non-universal 9 for a = 1, and non- 
power-law scaling implying an effective 9 — — oo for a > 1. 



which can be seen via saddle-point approximation. Tak- 
ing the normalization of the spins |Sj| — yfm into account 
and using Eq. (74), the main diagonal of C needs to be 
such that 



Cu — 1. 



(76) 



Then Eqs. (73), (75), and (76) can be solved self- 
consistently for the values of the N variables Hi. Eq. (76) 
is the saddle-point equation for the integrals in (71) and 
is valid as to — > oo at fixed N. 

At T = these equations are no longer well defined. 
Instead, as described above, it is a necessary condition 
for the system to be in a ground state for each spin to be 
aligned with its local molecular field H^, cf. Eq. (55). Due 
to the lack of metastability, for to — > oo this condition is 
also sufficient. As was first noted in Ref. 50, for T — > 
the variables Hi of Eq. (73) are equivalent to the rescaled 
amplitudes \Hi\/y/m of the local fields. 

It was suggested in Ref. 10 that the density p(A) of 
(real) eigenvalues of % _1 or A was useful in discussions 
of the physics of random spin systems. For example, 
it was shown that the smallest eigenvalue would van- 
ish at the critical temperature for N — > oo. For the 
case of the m = oo SK model in the thermodynamic 
limit, p(X) follows a Wigner semicircle. Decreasing the 
temperature from Tgc, where the first eigenvalue van- 
ishes, a fraction too ~ TV 2 / 5 of eigenvalues becomes zero 
as T — >• 0. 6 This corresponds to the contraction of the 
spin orientations into an mo-dimensional subspace in the 
limit of zero temperature, an effect reminiscent of the 
Bose-Einstein condensation in atomic systems. The be- 
havior of the eigenvalue density on cooling the system 
from high temperatures is illustrated for our Id system 
in the infinite-range regime a < 1/2 in Fig. 14. 



The spin integrations can be performed using integral 
representations of the S functions to yield 



ll~ ~ exp ' 



4tt 



E 



j3H l + lndet 



X 



where the susceptibility matrix x is defined by 



Xij — {-A )ij 
A-ij HiSij 



(71) 

(72) 
(73) 



The Hi in Eq. (71) are initially just new integration vari- 
ables introduced by the Fourier representation of the S 
constraints, but it will turn out that they have a profound 
physical meaning. The matrix A is obviously symmetric, 
since Jij — Jji. The correlation matrix, 



is identical to 



10 



C i 



c. 



1 



m 



T(A-%, 



(74) 



(75) 



B. Order parameter and spin-glass susceptibility 

The saddle-point equations for systems on hypercubic 
lattices in two and three dimensions (as well as in the 
SK limit) have been previously examined by Lee, Dhar 
and Young in Ref. 15. They considered the correlation 
function Cy of Eq. (75) and determined the Edwards- 
Anderson order parameter as 



3ea 



lim [C. 



(77) 



or as q^ A — [C?-] av , i ^ j for the SK model. Taking into 
account the scaling mo ~ of the number of non-zero 
spin components, they concluded that gj| A ~ in the 
ground state, i.e., that the order parameter vanishes in 
the thermodynamic limit. Similarly, defining 



Xs'g 



i,3 



(78) 



they inferred algebraically decaying correlations 



\Cl 



-dfi 



in the m = oo model, i.e., merely quasi 
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FIG. 14. (Color online). The eigenvalue density p(X) of the matrix A, see Eq. (73), for the Id m = oo model at a = 0.1. At 
high temperatures, the distribution follows Wigner's semicircle law 6 . Lowering the temperature shifts the whole distribution 
to the left, until for T — > mo eigenvalues {A a } vanish, while the majority (N — mo) of eigenvalues {At,} stay finite. In the 
thermodynamic limit the Wigner semicircle law is restored again also there, since mo/N — > 0. 



long-range order. This is in contrast to the findings of 
Ref. 7 for the case of the N — > oo limit being taken 
before the m — > oo limit. 

We believe, however, that one needs to consider the 
connected correlation function and on-site correlations to 
determine xsg and 5ea, and we will see that this leads 
to different conclusions. The basic idea is to separate 
contributions from the zero and non-zero modes. Con- 
sider Eq. (74) and factor out the cumulant part of 
the correlation function, 



Cij — 



1 



;s l ).(s 3 ) + l((s 



s,)- 



(S.MS,)). (79) 



Due to the symmetry of A and Eq. (75) the matrix C is 
also real and symmetric, and hence invertible. For the 
eigenvalues A of A and ui of C, we have the relation 



where a = 1, . . . , 
vanish as T — > 0, 
to the remaining 
Here, X x denotes 
component of the 
v x of A. 

Working at m 
tion that the a 1 
the form 



mo labels the mo eigenvalues A Q that 
and Af,, b = N — mo + 1, • • • , N refers 
N — niQ eigenvalues which stay finite. 



,l.h 



the x eigenvalue and v% is the k 
corresponding normalized eigenvector 

= mo and using the reasonable assump- 
component of the spin £ K™° has 



(83) 



the cumulant part of the correlation function can be iden- 
tified with 



(84) 



T 
A' 



(80) 



Using the spectral theorem, C has an orthonormal basis 
of real independent eigenvectors {u' l |n = 1, . . . ,N} for 
the set of eigenvalues {ui n g ]R}, which do not need to be 
distinct. 

The spectral decomposition reads 



N 



(if 



(81) 



where we use the outer product and the transpose de- 
noted by (-) T . Thus, with (v' l ) T ■ v 711 — S mn and for 
T < Tsq, the correlation matrix C has the eigenvalue 
decomposition 



being a function of the non- vanishing eigenvalues. It is 
natural to define the spin-glass susceptibility in terms of 
this connected correlation function, 



XSG 



Nm 



= — V r 2 =—X^ — 

ij 



(85) 



N ^ Xf 

b b 



where the Ab are the non-vanishing eigenvalues in the 
limit T — y 0. Furthermore, the Edwards-Anderson order 
parameter is then given by 
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FIG. 15. (Color online). Density plots for the parameters Hi in Eq. (73) for a = 0.1 (left) and a = 1.0 (right). In the limit 
T — > they correspond to the (rescaled) local fields \H.i\/i/m. As is seen here, Eq. (94) yields suitable starting values for the 
iteration at high temperatures, cf. the bold ticks on the Hi axis. For increasing values of a the distribution broadens. The 
normalization of the histograms is arbitrary. The system size used for this plot was TV = 1448. 



where the A a are the vanishing eigenvalues as T — > 0. 

The chosen formalism leads to a non-vanishing order 
parameter at T = 0. To see this, consider 



= 1 



(87) 



resulting from the normalization condition Cu = 1, 
Eq. (76). The finite eigenvalues \b scale to a constant 
as N — > oo and T — > 0, hence the second sum in Eq. (87) 
is roughly proportional to N^ 1 N 1 ^ fJ - = N - ^ and thus 
vanishes in the thermodynamic limit. Hence, q^A = 1 in 
this limit. The "zero" eigenvalues A a vanish as T — >• 
and as N — > oo. Assuming them to be proportional to 
T/NP, we conclude from Eq. (87) that p=l-fi. 



so that after multiplying with A 1 from the right and 
solving for the needed derivative we arrive at 



dA- 



dH, 



-A 



. x OA 



A~\ 



By virtue of Eq. (73) it is 
/ dA 

such that we finally find 
(dA- 1 



(90) 



(91) 



\9Hj 



EE( rl ),„M™»( rl ) m , 

n rn 

( A ~ 1 )ii( A ~ 1 )n = -( A ~ 1 )l- (92) 



C. Numerical approach 

We solve Eqs. (73), (75) and (76) iteratively using the 
Newton-Raphson method for systems with a < 0.85, see 
the discussion in Ref. 15. For larger a, this approach has 
some numerical instabilities leading to singular matrices 
in the course of the LU decomposition. We switched to a 
method using a QR decomposition similarly to the way 
Broyden's method 31 is usually implemented. However, 
in both cases we utilized the exact Jacobian to speed up 
calculations. Next we will show that this extra speed-up 
comes for free. According to Eqs. (73), (75) and (76) we 
can proceed by introducing the functions 

f i ({H k })=T(A- 1 ) ii -l (88) 

and solving for their zeroes. Taking into account A~ 1 A = 
1 we have 

51 dA- 1 . . , dA 

A + A~ —— = 0, (89) 



TABLE III. 


Realizations used for the T > calculations. 
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FIG. 16. (Color online). The evolution of eigenvalues A of the 
matrix A of Eq. (73) with temperature for a L = 1448 system 
with a = 0.1. The non-zero eigenvalues (black, dot-dashed) 
stay finite in the limit T — > (biggest are left away). The 
others (blue, solid) scale to zero as ~ T. Note that both axes 
have logarithmic scales. 

Hence the desired Jacobian reads 

( J fh = §§. = -T (A-% . (93) 

We start the iterations at high temperature where we 
expect 6 

Hi = T + 1/T (94) 

for the SK model, resulting in reasonable initial values 
also for the Id long-range system with a > considered 
here. This is illustrated in Fig. 15 for a = 0.1 and a = 1. 
For all systems considered here, we found a starting tem- 
perature of T\ = 10 to be a suitable and sufficient choice. 
The set of temperature points may be chosen in a geomet- 
rical fashion, such that the inverse temperature /3 = 1/T 
is distributed equidistantly. In all our simulations the 
lowest temperature was chosen to be Tf = 0.01. On 
decreasing the temperature, one may use the converged 
result of the previous calculation as a starting point. Al- 
ternatively, decreasing the temperature from to Tk+i 
(k > 0, Tk+\ < Tf.) a new guess for the values of the Hi 
can be obtained by using the differential equation 6 

^jf =-£(*"%. (95) 

3 

B l3 = {pC l3 f. (96) 

To perform the analysis outlined above in Sec. VB, we 
need to separate zero from non-zero eigenvalues. For fi- 
nite temperatures and finite systems, however, no eigen- 
values are exactly zero. Instead, there is a difference 



1 2 




0.0 0.2 0.4 0.6 0.8 1.0 



(7 

FIG. 17. (Color online). Spin-dimension exponent /i cal- 
culated from the eigenvalue density at finite temperatures 
(T > 0) as compared to the result reported in Sec. IV A from 
the ground-state computations (T = 0). For a = the SK 
model was considered directly. For the fully connected model, 
the results of the calculations at T = and for T > are well 
compatible. 



in scaling behavior between the "zero" eigenvalues that 
vanish proportional to T/N 1 ^^ and the other eigenvalues 
that scale to a constant. This is illustrated in Fig. 16. 
To automatically distinguish between the two types of 
eigenvalues, one might count all those as zero that fall 
below a chosen threshold at the lowest considered tem- 
perature. It turned out to be more reliable, however, to 
base the distinction criterion directly on the temperature 
scaling oc T of the "zero" eigenvalues. Determining the 
slope of logA(T) at the lowest considered temperatures, 
we counted those eigenvalues as scaling to zero whose 
slope was above 0.5. Even with this rather reliable cri- 
terion, however, there will always be a certain ambiguity 
as the slopes change quite continuously over the differ- 
ent eigenvalues, and a number of borderline cases always 
exists, cf. the example in Fig. 16. We do not find any 
signs of the number of zero eigenvalues changing with 
temperature. Instead, our results are compatible with 
all relevant eigenvalues starting to scale to zero as soon 
as T < T SG . 

The analysis of the eigenvalue density allows for an al- 
ternative method of calculating the spin dimension scal- 
ing exponent fi already discussed in Sec. IV A. To im- 
plement it, we determined the rank of the matrix M 
composed of the ground-state spin vectors. Instead, we 
could have extracted the local field values |Hj|, fed them 
into Eq. (73) and determined the number of zero eigen- 
values. These approaches are equivalent since the rows 
of M correspond to the null eigenvectors of A and the 
row and column ranks of a matrix are identical 6 . This 
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FIG. 18. (Color online). The eigenvalue distribution of the matrix A, cf. Eq. (73), at temperature T = 0.01 for different values 
of the power-law exponent a (N — 1448). Wigner's semicircle form holds for the infinite-range region < a < 1/2 in the 
thermodynamic limit. The red, bold line shows fits of the functional form (99) to the data. A positive exponent x results for 
a < 3/4. 



would not have allowed us to consider sufficiently large 
systems, however, since there an N x N matrix must be 
inverted, which is in contrast to the T = ground-state 
calculations where it was sufficient to determine the rank 
of the auxiliary, but smaller, to x N ground-state spin 
matrix M. The results for too extracted from the finite- 
temperature calculation are summarized in Fig. 17. The 
outcomes are mostly compatible with those of the zero- 
temperature approach for the case of the fully connected 
model. In contrast, the T — results for the diluted 
model systematically deviate from those of the fully con- 
nected model for a > 1 as discussed above in Sec. IV B, 
and some signs of this non-universality are already seen 
for a > 0.8. 



D. Distribution of eigenvalues 

Braun and Aspelmeier 58 suggested using the eigen- 
value spectrum for a more general understanding of scal- 
ing corrections for the case of the two competing limits 
N — » oo and m — > oo. Their analysis is valid for the sys- 
tem on a Bethe lattice, but some results might generalize 
to the model considered here. They discuss the ground- 
state energy e(m,N) = E/Nm per spin and spin com- 
ponent, which is argued to have two contributions: the 
ground-state energy in the limit N — > oo with m = to 
large and fixed, + \rriQ V + O(toq), and the additional 
energy required for forcing the N spins into an mo dimen- 
sional subspace. This second contribution is proportional 
to the required shift of the eigenvalue spectrum p(A) to 
push too eigenvalues to zero. Assuming the density at 
small A to scale as p(X) ~ X x , 10 the first too eigenval- 
ues occupy the interval [0, (mo/N) 1 ^ 1+x '>], such that the 
total energy is 
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FIG. 19. (Color online). The density exponent x as deter- 
mined by fitting a power-law, Eq. (99) , to the eigenvalue den- 
sities shown in Fig. 18. 



Minimizing with respect to the number too of spin com- 
ponents yields the scaling relation 



M = 



1 

y(x + 1) + 1 ' 



(98) 



For the SK model with x = 1/2 (the Wigner semicircle 
law) and y = l, 8 we therefore arrive at the observed 
/i = 2/5 as desired. This scaling should hold independent 
of lattice structure. 

We determined the exponent x of the density of eigen- 
values as A — > from fits of the functional form 



p(X) = <z(A + AX) X 



(99) 



e(m, N) 



eoo + cito v 



/TO0_\ 

V N ) 



i/(i+x) to the data. Here, the shift AA is required to take the zero 

• (97) 

eigenvalues into account. To perform the fits, eigenvalues 
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agree for a < 1/2 where b — z — 2/5. For larger values of 
a, however, b is consistently smaller than z. Therefore, if 
the corrections predicted here are present, they are sub- 
leading and cannot be resolved by our numerical analysis. 



E. Critical behavior 

We now turn to studying the behavior of the physical 
observables extracted from the solution to the saddle- 
point equations in the vicinity of the critical point. The 
analysis of the spin-glass correlation length, the Edwards- 
Anderson order parameter and the spin-glass susceptibil- 
ity allows us to compare our simulations to the theoretical 
predictions outlined in Sec. HID. 



Correlation length 



FIG. 20. (Color online). A comparison of the exponent z = 
(1 + /i)/(l +x) determined from the T > calculations (black 
diamonds) with the ground-state energy correction exponent 
b as determined in Sec. IV A (green circles). 



from a large number of disorder samples were accumu- 
lated, and the limit A — > was modeled by successively 
omitting more of the larger eigenvalues while monitoring 
the resulting estimate of x as well as the goodness-of-fit. 
Statistical errors on the ht results were determined us- 
ing a sophisticated jackknifing analysis 74 . Some example 
results are collected in Fig. 18. Fits of this functional 
form are possible for a < 3/4, where the vanishing of the 
phase transition is signaled by x = 0. This is expected 
since at the upper critical value a u — 3/4 the vanishing 
of eigenvalues at any finite temperature ceases to exist 10 . 
Collecting the results for all values of a considered, we 
arrive at the data shown in Fig. 19 which confirms our 
expectations of x — 1/2 for a < 1/2 and x — for 
a = 3/4. In view of the above expectations which are 
in-line with our numerical results, it is tempting to spec- 
ulate that x(a) = 3/2 — 2er, but we have not been able 
to substantiate this claim with a theoretical argument. 

Additionally, the authors of Ref. 58 suggest the scaling 
ansatz 



e(m, N) — e c 



(100) 



where F(x) is a scaling function. In the relevant limit 
of m — > oo before N — > oo, this implies a scaling of the 
ground-state energy according to 



e(m = oo, N) 



N~ 



(101) 



with z = \iy = (1 — /i)/(l + x), where we used Eq. (98). 
In Fig. 20 we show the correction-to-scaling exponent b 
for the ground-state energies as determined in Sec. IV A 
in comparison to z = (1 — fi)/(l + x) as determined from 
the results of fi and x for our system. Both exponents 



The spin-glass correlation function can be calculated 
from the spin-spin correlation function (75) via 



L 



E 

r \ ■ —r 

13 



(102) 

Note that here we use the algebraic graph distance = 
min(|i — j\,L — \i — j\) irrespective of whether the ring 
or chain geometry is considered. To arrive at the usual 
second-moment definition of the correlation length, we 
use the Fourier decomposition, 
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and plug it into Eq. (32), 

X° SG (0) 
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-, 1/(2,7-1) 



XSG^n 



(103) 



(104) 



Here, fc m ; n = 2n/L. In practice, we determine Gs G (r) 
per disorder realization from the saddle-point equations. 
For space efficiency, storing Gsc(r) is then preferable 
over storing Cy directly. Note that since we are using 
the disconnected correlation function here, the estimators 
(104) only represent the correlation length above Tsg- 
Close to criticality, we expect the scaling form 



L^^XitL 1 / 4 ), 1/2 < o < 5/8, 
L X{tl}l v ), a > 



(105) 



In the ordered phase, on the other hand, £1, diverges even 
more strongly with the system size 75 . As a consequence, 
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FIG. 21. (Color online). The spin-glass correlation length £l/L (for a > 5/8) and £/L"' 4 (for 1/2 < cr < 5/8), respectively, for 
different lattice sizes L. For a < 3/4 lines show a clear crossing point. The crossing vanishes at a — 3/4, where the lines for 
lattice sizes L = N > 362 lie approximately on top of each other below a certain temperature. 



UL/> 





L 


o 


512 


A 


724 




1024 


X 


1448 



T SG = 0.8(2) 
a= 1/4 
b = 0.08(5) 



cr = 0.54 Af- 



-10 12 
(T - T SG )L a 



o 
id 




T SG = 0.00(4 
a = 0.000(6 

b = 1 









L 


O 








i — 1 




o 


512 






A 


724 


lO 




+ 


1024 


O 




X 


1448 



(7 = 0.75 t% - 
I I I Jl 



0.0 0.1 0.2 0.3 0.4 0.5 

(T - T SG )L a 




0.1 0.3 

(T-T SG )L a 



FIG. 22. (Color online). Scaling collapse of the correlation length ratios £l/L and (,l/L v ^ a , respectively. For optimizing 
parameters, a general form £l = L b X(tL a ) was assumed. For a > 3/4 the transition temperature was fixed at Tsg = 0, since 
the according collapses yielded values around zero. The resulting parameters are summarized in Fig. 23. 



the curves for £ L /L (cr > 5/8) and £ L /L V / A (1/2 < cr < 
5/8), respectively, will cross in the vicinity of the critical 
temperature. 

Plots of £,l/L and $,l/L u / 4 resulting from the finite- 
temperature calculations with parameters summarized in 
Tab. Ill for three examples of cr are shown in Fig. 21. Be- 
low a u — 3/4, we find a crossing of the curves as shown 
for cr = 0.54. For a > 3/4, on the other hand, the curves 
only merge in the limit T — > 0. At the critical cr, we 
see a merging of the curves with an onset temperature 
scaling to as L —> oo. We note that the scaling of 
£,l/L"/ 4 is hard to observe numerically for a close to 1/2, 
where v — l/(2o p — 1) — > oo, such that a crossing point 
of £,l/L v / 4 is not visible for the system sizes considered 
here for a = 0.51. It is possible to extract estimates 
for the spin-glass temperature Tsg and the correlation 
length exponent v by re-scaling the data for different sys- 
tem sizes such that they collapse on the scaling function 
X(x). We used two complementary approaches for per- 



forming this collapse: method (a) consists of a joint fit of 
all data sets to a third-order polynomial approximating 
the scaling function in the chosen regime; method (b) is 
the collapsing procedure suggested in Ref. 76 which, in 
turn, is based on Ref. 77. In both cases, we performed 
the collapse on the logarithm of the actual data. This 
turns out to be necessary since, in particular for small 
°") spans many orders of magnitude in the range 

of temperatures considered here. In some cases, we also 
employed weights of the data points involved that decay 
exponentially away from the adaptively chosen value of 
the critical temperature. Statistical errors on the col- 
lapse parameters have been determined by a bootstrap 
sampling over the whole collapsing procedure 74 . In the 
region cr > 5/8, v was determined from the scale tL l l v 
of the abscissa. On the contrary, for the mean-field re- 
gion 1/2 < cr < 5/8, it was determined from the scaling 
^l/L"/ 4 of the ordinate. For the latter collapses, we find 
that the expected scaling of iL 1 / 4 of the argument of the 
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FIG. 23. (Color online). Spin-glass transition temperature Tsg and correlation length exponent v as estimated from collapses 
of the correlation length data according to the functional form Eq. (105). Left panel: transition temperature Tsg as compared 
to the prediction Tsg = \/3 — 4a of Eq. (51). Right panel: 1/v and the relations 1/v = 2a — 1 valid in the mean-field region 
1/2 < a < 5/8, cf. Eq. (34), and 1/v w (3 - 4<r)/2, cf. Eq. (52). For a > 3/4, where T SG = 0, we additionally show 1/v = -6 
with our estimates of 9 resulting from the defect-energy calculations, as well as the conjectured 1/v — —9 — a — 3/4. 




FIG. 24. (Color online). Scaling of the Edwards-Anderson order parameter ?ea- The left panel shows the dependence of qea 
on temperature for a — 0.1. For L — > 00, we expect the form (107) which is indicated by the solid blue line. Employing fits of 
the functional form (108) for a temperature interval T; < T' < Tf around the expected critical temperature (see inset of left 
panel), we determine Tsg and /3/v from the point where such fits work best, which is monitored by the quality-of-fit parameter 
Q, shown in the right panel. For a summary of results for all a see Fig. 25. 



scaling function is not very well reproduced, and we al- 
low for this exponent to fluctuate to accommodate scaling 
corrections. 

As illustrated in Fig. 22, these collapses work rather 
well over the whole range of a. The resulting estimates 
of the correlation length exponent v and the critical tem- 
perature are summarized in Fig. 23. The transition tem- 
perature is consistent with Tsg = for a > 3/4 and 
approaches Tsg = 1 as cr — > 1/2. In between, it is com- 



patible with the estimate Tsg ~ \/3 — 4a obtained in 
Sec. HID. As mentioned above, in the mean-field regime 
with a — ¥ 1/2+, finite-size corrections become very pro- 
nounced due to the divergent exponent v. This leads 
to rather strong fluctuations of v and Tsg as estimated 
from the collapsing procedures, cf. Fig. 23. In the right 
panel of this figure we also compare our result for 1/v 
extracted from collapses for a > 3/4 with —9 from the 
defect-energy calculations. In general, we find acceptable 
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agreement between zero- and finite-temperature calcula- 
tions. The observed systematic deviations give an indi- 
cation of the level of unresolved finite-size corrections. 
As a ~ > 1—, results for the diluted system start to sys- 
tematically deviate from those for the fully connected 
system due to the observed non-universality discussed in 
Sec. IV B. 

Comparing the estimates for l/v and T SG for the 
ring, line and resummed line geometries introduced in 
Sec. II B, we find complete consistency, cf. the example 
results for a — 3/4 collected in Table IV. 



2. Edwards- Anderson order parameter 

According to the discussion of FSS in our model, we 
expect 



/ L-^QitL 1 / 4 ), a<5/8, 

QEA ~ | L -P/»Q(t L l/^ u > 5/ 8 . 



(106) 



for temperatures in the scaling window. In the thermo- 
dynamic limit, for a < 1/2, f3 = l, 6 ' 7 while for a > 5/8, 
(3 is expected to remain close to unity, so that 



c/ea 



l-T/Tg G) T<T SG , 
0, T>T SG . 



(107) 



As is illustrated with the unsealed data in the left panel 
of Fig. 24, these expectations are borne out well by our 
results. In particular, the thus defined order parameter 
becomes unity as T — > 0, in contrast to the differently 
defined q^ A of Eq. (77) and Ref. 15. 

Right at Tsg, 9ea scales to zero. To extract Tsg and 
determine f3/is, we again employed scaling collapses. Due 
to the observed instability of the collapsing procedure, 
we also developed an independent approach based on the 
quality of power-law scaling. Since scaling proportional 
to 1_PI V is only expected at criticality, the critical point 
might be determined under the assumption that it co- 
incides with the temperature where power-law scaling is 
best observed. We hence performed fits according to the 
form 



5ea(T = T') = elf'", 



(108) 



for an interval of temperatures T < T' < Tf around the 
expected value of Tsg- If power-law scaling only occurs 
at T = TgG asymptotically, the quality-of-fit parameter 
Q should be maximized at this point, such that the infor- 
mation of both the critical temperature and the exponent 
fi/v can be extracted by this procedure. An example for 
this approach for a = 0.1 is shown in the right panel of 
Fig. 24. 

The overall results for the transition temperature Tsg 
and the critical exponent j5/v resulting from this anal- 
ysis are summarized in Fig. 25, together with the cor- 
responding results of a collapsing procedure. The esti- 
mates of the spin-glass temperature are consistent with 



Tsg = \/3 — 4er in the relevant regime and become con- 
stant at Tsg = 1 f° r G < 1/2, while they vanish for 
a > 3/4, as expected. From the present analysis, Tsg can 
be resolved with more precision than from the correlation 
length particularly in the mean-field regime a < 5/8. 
The exponent f3/v is consistent with the expectations 
summarized in Sec. HID, i.e., (3/v = 1/4 for a < 5/8 
and (3/v = (3 - 4cr)/2 for 5/8 < a < 3/4. The statistical 
precision of our determination, however, is not sufficient 
to rule out possible different scenarios and, in particular, 
to decide whether /3/v = (3 — 4cr)/2 might be exact in the 
non-mean-field regime. Again, statistical errors are cal- 
culated by an elaborate jackknifing procedure. As shown 
in Table IV, no significant deviations between the results 
for the different models of a Id geometry are observed. 



3. Spin-glass susceptibility 

We finally analyzed the scaling behavior of the spin- 
glass susceptibility as defined from the connected cor- 
relation function in Eq. (85). From the discussion in 
Sec. Ill D we expect scaling according to 



XSG 



L^CijtL 1 ^) a < 5/8, 
Uil v C(tL x l v ) cr > 5/8. 



(109) 



In contrast to the scaling of </ea it is possible here with- 
out reference to numerical derivatives to define a series 
of pseudo-critical temperatures from the locations of the 
maxima of the susceptibility, 

T^ ax) =T SG + cL- 1/l/ , (110) 

while the values of xsg at the maxima should then follow 

(111) 



Fits of the corresponding forms to the data for cr = 0.1 
are shown in Fig. 26. We find, however, that the resulting 
parameter estimates are afflicted by very strong finite- 
size corrections. In particular, the resulting estimates of 
1/v are far off from our theoretical predictions as well as 
the results from the analysis of the correlation length. 

The presence of strong corrections in the scaling of xsg 
is well-known from studies, e.g., of the Ising spin glass. It 



TABLE IV. Results for different choices of the geometry of 
the model, which were introduced in Sec. II B. There are no 
significant differences for the value a — 3/4 checked here. 



quantity 


ring 


line 


summed line 


A 4 


0.372(6) 


0.374(8) 


0.371(7) 


1/v 


0.000(5) 


-0.01(1) 


-0.03(5) 


Tsg 


0.00(4) 


0.04(6) 


0.01(16) 


<?EA 










0.01(4) 


-0.02(5) 


0.00(2) 


Tsg 


0.06(19) 


-0.09(27) 


0.00(13) 
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FIG. 25. (Color online). Transition temperature Tsg and critical exponent fi/v as extracted from the Q-maximization procedure 
described in the main text and in Fig. 24. As a comparison we extracted these quantities from a data collapse as well. 
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FIG. 26. (Color online). Scaling of maxima of the spin-glass 
susceptibility xsg in a showcase example with a = 0.1. 



has been suggested 78 that modified scaling forms incor- 
porating scaling corrections might contribute towards re- 
solving such corrections and the proposed extended scal- 
ing forms have been successfully applied to the Ising spin 
glass 79 . In particular, one problem of the scaling form 
(109) is that it cannot reproduce the observed behav- 
ior Xsg — > 1 as T — > oo, cf. Fig. 26: assuming that 
C(x) ~ x a for x 3> 1, asymptotic size independence of 
the data at large T requires that a = —7 and hence 
LPtl v C(tL 1 l v ) -s- as T -> 00. While this is not in con- 
tradiction to scaling theory as the assumed scaling form 
should only apply in the critical region, having a scaling 
form consistent with the behavior as T — ¥ or T — > 00 
might allow to extend the scaling regime or, equivalently, 
reduce the observed finite-size corrections. A modified 
scaling form that serves this purpose is given by 78 

Xsg = (LTyl u C[(LT) l l v t], (112) 

which is compatible with \sg — >■ 1 as T — s- 00. We 



used this extended scaling form to perform collapses of 
the finite-size data for the spin-glass susceptibility. Even 
though some scaling corrections are implicitly included 
in Eq. (112), these collapses are found to be rather un- 
stable and, hence, sensitive to the choice of starting val- 
ues for the parameters and the range of data points to 
be included for each lattice size. In view of these un- 
certainties, we found it impossible to extract all three 
parameters, T$q, 7/1/ and 1/v reliably from a single col- 
lapsing procedure. We hence decided to keep Tsg fixed 
at the theoretical prediction TgG = \/3 — 4a which, as 
is shown in the results of Fig. 24 for the order param- 
eter, is well compatible with our numerical results. An 
example collapse is shown in the left panel of Fig. 27. 
The right panel shows our resulting estimates of 7/z/ for 
a < 0.8. These are roughly compatible with our expec- 
tations of 7/1/ = 0.25 for a < 5/8 and ajv = 2a — 1 
for (7 > 5/8. For a > 3/4, we do not find stable col- 
lapses with reasonable parameters which we attribute to 
the fact that, there, Tsg = 0, but our data only reach 
down to T m j n = 0.01. The resulting values of 1/v are 
strongly fluctuating and hence not useful as reliable esti- 
mates of this quantity. An alternative collapsing exercise 
using a plot of xsg as a function of £/L, which should 
have the theoretical advantage of involving only a single 
adjustable parameter j/v did, unfortunately, not lead to 
more reliable results. 



VI. CONCLUSIONS 

We have used extensive numerical simulations together 
with a number of phcnomcnological scaling arguments 
to give a rather comprehensive account of the critical 
behavior of the one-dimensional spin-glass model with 
power-law interactions in the limit of an infinite number 
of spin components m. Compared to the more familiar 
case with m < 00, we find a number of remarkable differ- 
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FIG. 27. (Color online). Top: A collapse of xsg accord- 
ing to the extended scaling form (112) for a = 0.1 with 
adaptively determined parameters j/u = 0.24 and v — 0.37 
(Tsg = 1). Bottom: estimates of 7/1/ resulting from an adap- 
tive collapsing routine together with the theoretical expecta- 
tions f/v = 0.25 for a < 5/8 and a/v — 2a -— 1 for a > 5/8, 
respectively. 



ences which are, in part, related to the shift in the lower 
and upper critical dimensions of the model. 

The lack of metastability in the model allows to per- 
form quasi-exact ground-state calculations. The result- 
ing defect energies are well described by a long-range 
stiffness exponent 9 = 3/4 — a. This relation can also 
be deduced from scaling arguments, but we have not 
been able to provide a more rigorous derivation. This 
relation results in an upper critical <r u = 3/4, where 
finite-temperature spin-glass transitions first disappear. 
On lowering a, mean- field behavior sets in at <r; = 5/8. 
These critical interaction ranges arc different from the 
(7/ = 2/3 and a u = 1 found for spin glasses with finite 
771. 17 For hypercubic lattices it has been speculated that 
di = d u if 777, — > 00. 14,15 According to our analysis of the 
Id model theoretical as well as numerical evidence exists 



showing that the upper and lower critical a are well sep- 
arated. We also investigated the distribution of ground- 
state energies, and it is found to be Gaussian for the 
full range of a, again in contrast to the Ising case where 
non-Gaussian distributions where found in the mean-field 
regime 20 . Sample-to-samplc fluctuations are trivial with 
Of — 1/2 for a > 1/2, but cross over to a value consis- 
tent with Qf = 1/5 for a < 1/2 as conjectured for the 
SK model in Ref. 63. 

In Ref. 26 it has been suggested to study a diluted 
version of the Id long-range spin-glass model to reach 
even larger system sizes. As we point out here, however, 
the two models are not in the same universality class for 
a > 1, where the diluted model becomes equivalent to 
a short-range Id system. Right at <r = 1, critical expo- 
nents depend continuously on the average coordination 
number z. Also, samplc-to-sample fluctuations are triv- 
ial for the diluted model with Qf — 1/2 for all a, an 
effect anticipated for the Ising case 62 . Additionally, we 
observe more pronounced scaling corrections for the di- 
luted model even in the regime 1/2 < a < 1 such that 
we cannot find an advantage for numerical simulations 
in the larger system sizes reachable through the dilution. 
Additionally, we compared different realizations of the Id 
geometry using rings and chains with and without Ewald 
summation of interactions. As scaling corrections for this 
class of models are pronounced, one could hope that some 
variant of the model leads to a substantial reduction in 
corrections. For zero temperature (cf. Fig. 8) as well as 
for the critical behavior at T > 0, however, we find no 
significant differences in the FSS of the different model 
variants considered, cf. the data collected in Table IV. 

The ground-state calculations are complemented by re- 
sults from an iterative solution of the saddle-point equa- 
tions resulting in the m — > 00 limit, yielding access to the 
order parameter, spin-glass susceptibility and correlation 
length. In contrast to Ref. 15, we argue that using appro- 
priate definitions of the basic observables the model does 
show true long-range order in the low-temperature phase, 
even for the order of limits m — > 00 before N — > 00 nat- 
urally taken in numerical studies. The critical exponents 
v, j3 and 7 numerically determined from this approach 
are consistent with our theoretical arguments in the full 
range of a. The critical exponents in the non-mean- 
field regime have been hard to determine with precision. 
However, one of the surprises is the utility of the simple 
approximate RG scheme first suggested by McMillan 47 
which seems to work quite well over the entire non-mean 
field region, 5/8 < a < 3/4. 

As an aside, we find clear-cut evidence of the recently 
suggested exactness of mean-field theory for spin-glass 
models in the regime a < 1/2, 30 where we see non- 
universal properties such as the average ground-state en- 
ergy to be independent of the interaction range a. 

The m — » 00 model studied in this paper is an inter- 
esting model in its own right, partly because it is one 
of relatively few models known which have a failure of 
hyperscaling. The phenomenon of dimensional reduction 
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occurs in its short-range d-dimensional version which sug- 
gests that there might be some elegant supersymmetry 
in the model, but this has yet to be discovered. But our 
chief motivation in understanding this model was to clear 
the ground for our future 1/m expansion study of spin 
glasses. 
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